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Abstract. We present a study comparing the energy carried away by a 
coronal mass ejection (CME) and the radiative energy loss in associated flare 
plasma, with the decrease in magnetic free energy during a release in active re- 
gion NOAA 10930 on December 13, 2006 during the declining phase of the solar 
cycle 23. The ejected CME was fast and directed towards the Earth with a pro- 
jected speed of ^1780 km s" 1 and a de-projected speed of ~3060 km s _1 . We 
regard these as lower and upper limits for our calculations. It was accompanied 
by an X3.4 class flare in the active region. The CME carried (1.2-4.5) x 10 32 erg 
(projected-deprojected) of kinetic and gravitational potential energy. The es- 
timated radiative energy loss during the flare was found to be 9.04 xlO 30 erg. 
The sum of these energies was compared with the decrease in measured free 
magnetic energy during the flare/CME. The free energy is that above the mini- 
mum energy configuration and was estimated using the magnetic virial theorem. 
The estimated decrease in magnetic free energy is large, 3.11 xlO 32 erg after the 
flare/CME compared to the pre-flare energy. Given the range of possible ener- 
gies we estimate that 50-100% of the CME energy arose from the active region. 
The rest of the free magnetic energy was distributed among the radiative energy 
loss, particle acceleration, plasma and magnetic field reorientation. 
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1. Introduction 

Sunspots are the harbor for solar magnetic fields. Energetic phenomena like solar flares 
occurring in sunspots derive their energy from the active region magnetic fields. It is 
generally believed that coronal mass ejections (CME) which are associated with active 
regions obtain the majority, if not all of, their energy from these magnetic structures. The 
magnetic energy stored in the twisted magnetic fields is released within a few minutes 
during these energetic events and is of the order of 10 32 erg. The majority of the energy 
released is carried away by the CME mostly in the form of kinetic energy and partly 
as gravitational potential energy. The rest of the energy is mostly utilized in particle 
acceleration, plasma heating, radiation and magnetic field re-orientation (Webb et al. 
1980, Canfield et al. 1980, Emslie et al. 2004). 

Most of the energetic processes occur in the corona. To understand these processes 
it is essential to measure the magnetic energy and its temporal evolution. However, 
measurement of the coronal magnetic field is difficult because of observational constraints 
(Solanki et al. 2003, Lin et al. 2004). Currently, only the magnetic field measurements 
are accurately obtained at the photospheric levels with sufficient resolution and cadence. 
With the available magnetograms and using many techniques it is, however, possible to 
estimate the available magnetic energy in the corona with some uncertainties. 

There are many methods to estimate the available magnetic energy in the active re- 
gion corona. These are: (1) the minimum-current corona model (Longcope et al. 2007, 
Kazachenko et al. 2009); (2) the non-linear force free field (NLFFF) extrapolation and 
using the volume integral (Srivastava et al. 2009); (3) the Poynting flux estimation 
(Ravindra, Longcope and Abbett 2008); and (4) the magnetic virial theorem (Sakurai 
1987, Mctcalf et al. 1995, Mctcalf, Leka and Mickey, 2005, Venkatakrishnan and Ravin- 
dra, 2003). Among these, the energy estimation using the minimum-current corona model 
requires only the line-of-sight magnetic field measurements and it provides only the energy 
involved in the flare. All other techniques require the vector magnetic field measurements 
taken at the photosphere or chromosphere. 

The energy estimation using the Poynting flux requires not only the magnetic vec- 
tor field measurement, but also the magnetic footpoint velocity. The remaining two 
techniques rely on vector magnetograms obtained at one level, preferably in the chromo- 
sphere where the magnetic field is force-free. Up until very recently, the vector magnetic 
field measurements have been routinely made at the photospheric level. With the ad- 
vent of better computing facilities and readily available photospheric data, researchers 
have developed algorithms and applied them to photospheric vector magnetograms to 
mimic the chromospheric data (Wiegelmann, Inhester and Sakurai 2006). The procedure 
is known as the "pre-processing" of the vector field data. This technique was further 
improved when the chromospheric H Q data became available (Wiegelmann et al. 2008). 
The resulting magnetograms are used as lower boundary data for the NLFFF extrapola- 
tion to compute the coronal magnetic fields (Derosa et al. 2009). After computing the 
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coronal magnetic fields it is a relatively simple task to estimate the free energy available 
within the volume (Schrijver et al. 2008). 

The magnetic virial theorem can be applied to estimate the magnetic energy provided 
the field measurements are made at a level where the magnetic force and torque vanish 
(Molodenskii 1969, Sakurai 1987, Metcalf et al. 1995). In the past, many researchers have 
used the magnetic virial theorem to compute the free magnetic energy (Gary et al. 1987, 
Metcalf et al. 1995, Metcalf, Lcka and Mickey 2005). Metcalf et al. (1995) and Metcalf, 
Lcka and Mickey (2005) used the chromospheric vector magnetic field data observed in 
the Na I absorption line and these are considered to be close to the force-free condition 
(Metcalf et al. 1995). Using the chromospheric vector magnetograms they not only 
verified the force-free conditions, but they also provided a measurement of the magnetic 
energy. The magnetic free energy was measured in NOAA AR 10486 (a super active 
region) on 29 October 2003 before and after the flare (Metcalf, Leka and Mickey 2005). 
They found an increase in magnetic free energy after the flare compared to the pre-flarc 
energy level. 

Space weather prediction is a major challenge to the solar community as the major 
eruptive events on the Sun can produce large geomagnetic storms. Since the eruptive 
events derive their energy from the active regions, it is essential to measure the available 
free energy in the active region volume. It is also essential to validate this based on 
the energy carried away by the eruptive processes (Emslie et al. 2004). In this paper, 
we compare the drop in free magnetic energy with the energy content of the CME and 
thermal plasma during the X3.4 class flare produced by AR 10930. 



1.1 Active Region 10930 

The active region NOAA 10930 produced several X-class and M-class flares during its 
passage across the Sun in December 2006. It was well observed by Hinode and SOHO 
whose data sets are not influenced by seeing and Earth's diurnal effects. This active region 
has been extensively studied for the measurement of twist (Su, et al. 2009), rotation (Min 
and Chae 2009), penumbral filaments (Tan et al. 2009), flares (Isobe et al. 2007) and 
helicity (Magara and Tsuneta 2008). Space based observations, free from seeing induced 
polarization, better techniques in measuring the Stokes vectors, improved sensitivity of 
the instrument, good spatial resolution, improved inversion and ambiguity resolution 
techniques, and techniques to mimic the vector field measured at the photospheric level 
to closely match the vector field measured at the lower chromosphere, all provide good 
data sets for the measurement of free magnetic energy in AR 10930 before and after the 
X3.4 class flare occurred on Dec 13, between 02:14 and 02:57 UT. 

In this paper, we present the study of the temporal evolution of magnetic energy 
computed over a period of 5 days surrounding the flare and associated Earth-directed 
CME. We used the Hinode/SP vector magnetograms to estimate the magnetic free energy 
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available in the active region during this period. We applied a magnetic virial theorem 
algorithm to the pre-processed vector magnetic field measurements. In Section 2, we use 
vector magnetograms computed from Hinodc/SP to estimate the amount of accumulated 
magnetic free energy during the evolution of AR 10930 from December 9-14, 2006. This 
time frame includes the flare and CME occurred on December 13. To better understand 
the decrease in free energy from this active region, we estimate the amount of kinetic 
and gravitational potential energy carried by the CME using the time series of LASCO 
images. In Section 3, we present the evolution of magnetic free energy using the magnetic 
virial theorem over the time period. A comparison of the decrease in magnetic free energy 
after the event with the energy carried away by the CME and the radiative energy loss 
during the flare is performed. In the last Section, we conclude with discussions on the 
energy partition of the CME and flare. 

2. Observational Data and Analysis Methods 
2.1 Hinode/SOT 

The solar optical telescope (SOT) on board the Hinode satellite (Kosugi et al. 2007) 
makes spectro-polarimetric measurements at a resolution of 0.3" (Ichimoto et al. 2008). 
The active region maps have been produced by spatial scanning and the Stokes I, Q, U 
and V spectra are obtained in the Fe I 6301.5 A and 6302.5 A absorption lines. The 
spectro-polarimeter is capable of making the raster scan of the active region in the fast 
mode as well as in the normal mode. In the fast mode, spatial resolution along the 
slit direction is 0.295" and in the scanning direction is 0.317"/pixel. The fast mode 
scanning takes typically 1 hr and the normal mode takes 2 hr to complete the I, Q, U 
and V measurements over the active region. Our data set consists of only the fast mode 
data. The Stokes data were obtained from 9-14 December 2006. We chose the period at 
which the active region was located in the solar longitude range of E30 to W30. During 
this period 24 magnetograms were obtained from the fast mode. The Stokes signals 
were calibrated using the standard solar software pipeline for the spectro-polarimetry. 
In order to get the complete information on the vector magnetic fields, the obtained 
Stokes vectors have been inverted using the Milne-Eddington inversion (Skumanich and 
Litcs 1987, Lites and Skumanich 1990, Lites et al. 1993). The 180° ambiguity was 
resolved based on the minimum energy algorithm (Metcalf 1994) implemented by Leka, 
Barnes and Crouch (2009). The resulting magnetic field vectors have been transformed 
to heliographic coordinates as described by Venkatakrishnan, Hagyard and Hathaway 
(1988). The resulting B z (line-of-sight field) and transverse field strength has 1-cr error 
bars of 8 G and 30 G respectively. These vector field data have been used in computing 
the available magnetic energy. 

SOT/SP makes the vector magnetic field measurement at the photosphcric level where 
the force-free assumption is invalid (Metcalf et al. 1995). Since the measurements of the 
vector magnetic fields at the chromospheric/coronal level are not available routinely, we 
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have used the field measurements made at the photospheric level. Wiegelmann, Inhcstcr 
and Sakurai (2006) have developed a method of pre-processing to mimic the photospheric 
vector magnetic field data so that they appear to be force-free. We have applied the 
pre-processing method described in Metcalf ct al. (2008) to the vector magnetograms 
to minimize the effect of the magnetic force and torque in the data. The procedure 
minimizes the 2-D functional as in their Equation (7). This function consists of 4 terms 
and each constraint is weighted by a undetermined factor. The first and second term 
in their equation corresponds to force balance and the torque-free condition. The last 
two terms correspond to optimized boundary condition and smoothing parameters. The 
resulting pre-processed data set is close to force free and torque free (see Section 3.1) but 
not necessarily in each and every point in the magnetogram. 

In order to estimate the magnetic energy for the force-free media we have used the 
magnetic virial theorem (Chandrasckhar 1961, Molodenskii 1969, Aly 1984 and Low 
1985), given by 



where, E is the magnetic energy, B x and B y are components of the horizontal mag- 
netic fields and B z is the vertical component of the field. The free magnetic energy is the 
excess energy above the minimum energy state (energy content of the potential field). 
The free energy Ef ree can be estimated by taking the difference in the total available 
energy and the magnetic potential energy as 



where Exot is the total accumulated magnetic energy and Ep ot is the magnetic po- 
tential energy. While computing the energy using the virial theorem we opted for the 
center of FOV of the magnetogram as the origin. The potential magnetic fields have been 
derived from the B z component by using the Fourier transform method (Alissandrakis, 
1981). The routine fff .pro available in the solar software pipeline has been adopted here 
to compute the potential component of the magnetic field, B xp and B yp . This routine 
is written based on the constant a force-free field method of Alissandrakis (1981) and 
Gary (1989). The pre-processed B z component is used in the computation of B xp and 
B yp . Using B z , B xp , B yp and Equation (1) we have computed the energy of the mag- 
netic potential field. For the computation of total magnetic energy we have used the 
pre-processed vector magnetograms. 




(1) 



Efree = E T ot — E Pot 



(2) 
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Figure 1. A halo coronal mass ejection observed by the LASCO C2 coronagraph at 02:50 UT 
(left) and C3 at 03:18 UT (right) on December 13, 2006. These are running difference images 
where the previous image has been subtracted from the present. The white circle at the center 
of each image represents the solar surface and the solid grey circle represents the occulting disk. 
In C3 the so-called "snow storm", due to high-energy particles striking the LASCO detector, is 
visible. These later saturate the images in C2 and C3 to the point where no measurements are 
possible. 

2.2 SOHO/LASCO Data 

To analyze the associated coronal mass ejection we have utilized space-based corona- 
graphs. The Large Angle Spectroscopic Coronagraph (LASCO) on board SOHO (Brueck- 
ner et al. 1995) currently consists of two coronagraphs that block out the majority of 
the bright photospheric light from the Sun to reveal the surrounding corona. The C2 
coronagraph has a field of view of 2.0-6.0 R© and a cadence of approximately 30 min- 
utes, while C3 has a field of view of 3.7-30 R© and a cadence of ^50 minutes (where, 
R Q corresponds to 1-solar radius). Both are white light telescopes, i.e. they observe the 
Thomson scattered light from free electrons in the plasma comprising the corona. Den- 
sity enhancements in the corona such as coronal mass ejections (CMEs) can be identified 
with relative ease, and the distance from the Sun and mass of the CME can be estimated. 
LASCO has shown to be highly successful at identifying and tracking CMEs, with well 
over 10 4 detecte d to date (see Yashiro et al. 2004, Howard et al. 2008 or the CDAW 
CME catalog at |http://cdaw.gsfc.nasa.gov/C ME_list). 

By applying well-established assumptions about the observed CME we may identify 
its distance evolution and mass, from which an estimate of its kinetic and gravitational 
potential energy may be produced. Distance measurements are obtained by measuring 
the location of the leading edge of the CME relative to the Sun. Most workers choose a 
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fixed direction or single location on the CME front at which to make these measurements 
(e.g. Yashiro et al. 2004). Mass is estimated by measuring the intensity across the entire 
area enclosing the CME and converting it to a density using Thomson scattering theory 
(e.g. Billings 1966). 

Projection effects play an important role in the observed structure of the CME (e.g. 
Howard and Tappin 2008). LASCO images are projected into the plane of the sky, and 
so the closer the CME is to the Sun-Earth line, the more heavily projected its images will 
be. For mass calculations the angle to the plane of the sky is an integral component of the 
coefficients involved in the Thomson scattering physics (angle \ in Billings (1966)), and 
so the direction of propagation can be easily included in the calculations. For identifying 
3-D distance properties we may apply spherical geometry provided we have an idea of 
the direction of propagation of the CME. For the analysis in the present study we apply 
the technique of Howard et al. (2007, 2008) who used the location of the associated solar 
surface eruption as the direction of propagation. In this case, this is the flare in AR 
10930, which was located at 6°S, 23°W at the time of the eruption. It is important to 
note that the location of the flare is not an excellent indicator of CME direction as it has 
been well established that flares are more usually associated with the footpoint of a CME 
(e.g. Harrison 1986). In the absence of other 3-D directional information, however, the 
location of the flare is a reasonable approximator. 

Because of these and other uncertainties we consider two values for speed and mass 
and consequently for the energies: the projected values (unchanged from the original 
height-time measurements); and de-projected values following the analysis below. These 
two values may be regarded as lower and upper limits to the calculated values of speed, 
mass and energy, respectively. 

Following Howard et al. (2007) we use the following equation for 3-D distance: 



— = sin a cot e + cos a, (3) 
it 

where R is the radial distance of the measured point (on the CME) from the Sun in AU 
and a is the angle subtended by the measured point at the Sun, given by 

cos a — cos 9 cos 0, (4) 

where 9 and <j> are the latitude and longitude of the direction of propagation vector, 
which in this case is the location of the flare site. The angle e is the elongation, which is 
easily obtained via the application of two simple assumptions about CMEs observed by 
coronagraphs. These are the Point P approximation (e.g. Houminer and Hewish 1972) 
and that e is small. This enables the simple conversion r (AU) = e (rad), where r is the 
projected distance of the measured point (Howard et al. 2007). To determine the mass 
we use the standard codes for mass calculation available via the Solarsoft software suite 
(including eltheory.pro). These codes are based on the theory of Billings (1966) (page 
150) and include the \ term for direction. 
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Figure 2. A map of the vector magnetic field of the active region 10930. The transverse field 
vectors are overlaid upon the B z component of the magnetic field. The black (white) colors 
represent the south (north) polarity regions. 

The CME associated with the flare on December 13 was Earth-directed event and 
hrst observed by LASCO C2 at 02:54 UT on December 13. Assuming a constant speed 
its projected onset was 01:49-02:19 UT (de-projected-projected), within the 8 hours 
between the vector magnetic field measurements. Figure Q] shows LASCO C2 and C3 
images of the CME. Shortly after launch, energetic particles from the flare bombarded 
SOHO resulting in the saturation known as a "snow storm" . The bombardment intensified 
reaching a point where the LASCO cameras were almost completely saturated by around 
04:00 UT. As a result, we were only able to obtain mass measurements from single C2 
and C3 images (at 02:54 and 03:18 UT from Figures[T]respectively) and could only obtain 
height (or distance) - measurement until 04:18 UT. 

3. Results 

The active region NOAA 10930 appeared on the Eastern limb of the Sun on December 5, 
2006 in the southern hemisphere at a latitude of 5°, with two polarities that were almost 
perpendicularly aligned with the East- West direction. It was visible on the solar disk until 
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Figure 3. Scatter plots of .B^top left), By (top right) and B z (bottom) between the photospheric 
and pre-processed vector field data. 



December 17 when it crossed the West limb of the sun (see: http://www.solarmonitor.org/index.php). 
Figure [2] shows the sample image of the vector magnetogram with the transverse field 
vectors overlaid upon B z for the Dec 11 observation. The vertical magnetic field strength 
in the South (negative) and North (positive) polarity sunspot is in excess of -4000 and 
2300 G respectively. 

In the vector field map (Figure it is clear that the field lines are twisted. Both 
the polarities are located close to each other and the boundary between the two is highly 
sheared. The North polarity region is emerging in the Southern part of the existing South 
polarity region. As it emerges it rotates in the anti-clockwise direction and the transverse 
fields are twisted in the clockwise direction. 



3.1 Magnetic Energy 



The pre-processed vector magnetograms have been used to compute the magnetic energy. 
Because of the pre-processing, the 3-components of the magnetic fields have been mod- 
ified. We compare the pre-processed vector field data with the original data. Figure [3] 
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Figure 4. The total magnetic energy (*), the magnetic potential energy (o) and the magnetic 
free energy (A) available in the active region is plotted as a function of time. The vertical dashed 
line represents the onset of the X3.4 class flare. 



shows the pixel-by-pixel scatter plot of the same for 3 components of the vector magnetic 
fields. We have kept the B z component of the magnetic field close to the original as sug- 
gested by Metcalf et al. (2008). However, the pre-processing has changed the transverse 
fields significantly. In order to make sure that the pre-processed vector magnetograms 
show the force- freeness assumption, we computed the force balance parameter (e force) 
and torque balance parameter {(-torque) as suggested in Wiegelmann, Inhester and Saku- 
rai (2006). We found that the force and torque balance parameters are much smaller than 
unity. For example, the computed ef orce — 0.002 and ttorque = 0.037 for the December 
11, pre-processed magnetogram data observed between 03:10 and 04:19 UT. Similarly, 
the computed e/ orce = —0.003 and ttorque — 0.03 for the December 11, between 08:00 UT 
and 09:30 UT. This value is much smaller compared to the e/orce = 0.17 and ttorque 
= 0.37 for the original Hinode magnetograms obtained on December 11, between 03:10 
and 04:19 UT. A similar value is obtained for the magnetogram taken on December 11, 
between 08:00 UT and 09:30 UT. 

Figure 0] shows the plot of the total energy (*), magnetic potential energy (o) and 
free magnetic energy (A) computed using Equations (1) and (2) as functions of time. 
The uncertainty in measuring the total magnetic energy involves (1) those from the mea- 
surement in the field and (2) those due to the flux imbalance term in the virial theorem. 
The uncertainty in the magnetic field measurement has been propagated through the 
magnetic virial energy equation and we have estimated the uncertainty in total magnetic 
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energy. The term 1/87T f sv B 2 (x.n)cPx is neglected in the virial energy estimation (see 
Metcalf et al. 2008). This introduces an uncertainty in the total energy estimation. We 
estimated the uncertainty introduced by the neglected term in the virial energy by making 
the assumption that the net flux is spread over the hemisphere uniformly. We have taken 
the radius of the region of interest as x in the equation and the area as the surface area 
of the hemisphere. The final uncertainty is \J erroriX) 2 + error (2) 2 , where error (1) and 
error(2) are the uncertainties due to reasons (1) and (2) respectively. The uncertainties 
obtained from these two are small and shown in the plot. There is one more systematic 
uncertainty that occurs due to the fact that the field is not force-free in the photosphere. 
We have used the pre-processing technique to make it more closely approximating the 
force-free condition. However, there are still residual net forces. These residuals produce 
a systematic uncertainty in the virial energy estimation. Metcalf et al. (2008) model 
shows that the net Lorentz force underestimates the virial energy by about 10%. By 
taking this into account our estimated total virial energy goes up by 10% and hence 
also the free energy. This systematic underestimation shifts the total energy and free 
energy curve proportionately to each other. However this does not change the difference 
in energy between before and after the flare significantly. From the plot it is clear that 
the free magnetic energy increases with time. However, the increase in the free energy is 
large on December 11, at 23:10 UT to December 12, at 20:30 UT, and the free magnetic 
energy accumulated over a period of 20 hrs is estimated to be about 3.7xl0 32 erg. The 
last vector magnetogram before the X3.4 class flare was obtained on December 12, at 
20:30 UT. The next vector magnetogram was obtained on December 13, at 04:30 UT. So, 
there is a gap of about 8 hrs between the magnetograms. During this period we observe a 
decrease in the magnetic free energy of about 3.11 xlO 32 erg. This drop in free energy is 
clearly seen in the plot after the X3.4 class flare (shown as the dashed vertical line for the 
onset time of X3.4 class flare) has occurred. Afterwards the free energy again increases 
with time. 



3.2 Energy Carried Away by the Coronal Mass Ejection 

In order to compare the available magnetic free energy with the energy carried away by 
the CME, we have estimated the mechanical energy in the CME. The mass of the CME, 
averaged between the LASCO C2 and C3 images and assuming both the projected value 
and that at an angle of 67° from the sky plane, was (8.6-9.3) x 10 15 g. Its speed, obtained 
from a least squares linear fit through the distance-time plot and de-projected using the 
location of the flare and Equation ([3]), was 1775-3060 km s . This is a fast and massive 
CME but not physically unreasonable. CMEs have been known to obtain such speeds 
early in their evolution but are expected to decelerate rapidly through the heliosphere. 
Given that the travel time from the Sun to the ACE spacecraft (judging by the arrival 
of a forward shock at ACE on December 14 around 13:50 UT) was just under 1^ days, 
such a deceleration was likely. It is noteworthy that following a sudden commencement 
at 14:15 UT on December 14, a strong geomagnetic storm followed achieving a maximum 
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Figure 5. GOES-11 X-ray flux for 5 days starting from Dec 9, 2006. The upper and lower 
curves in the plot indicates the integrated full-Sun soft X-ray flux in the 1-8 and 0.5-4 A band 
passes respectively. The 1-8 A band pass curves give the X-ray flare index. 

Kp of 8+ and minimum Dst of -146 nT. Such activity is indicative of the arrival of a 
large fast CME at the Earth. 

The resulting kinetic energy of the CME was therefore (1.35-4.35) x 10 32 erg. We 
calculated the gravitational potential energy at each distance from the Sun but it never 
exceeded 4% of the kinetic energy. Including the gravitational potential energy we esti- 
mate the total energy to be (1.4-4.5) x 10 32 erg. This is (0.45-1.45) times the available 
magnetic free energy. 



3.3 Energy Content of Thermal Plasma 

Geostationary Operational Environmental Satellite (GOES) provides disk integrated X- 
ray flux in two wavelength bands, (1) 1-8 A and (2) 0.5-4 A (Kahler and Kreplin, 1991). 
The 1-8 A bandpass curves give the X-ray flare index. Figure [5] shows the GOES X-ray 
flux obtained from GOES-11 for the two channels over a period of 7-days around the event 
of interest. The data from the two channels were used to estimate the radiative soft-X-ray 
energy. We first subtracted the background flux taken on December 10, 2006 between 
19:09 UT and 23:12 UT. Using the solar software pipeline, the temperature of the plasma 
has been estimated by the taking the ratio of the two channels (Thomas et al. 1985). 
Using the synthetic solar spectrum, the computed temperature can be converted into an 
emission measure. By integrating this over the flare duration we obtained the radiated 
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energy by the thermal plasma (integrated between 02:14 UT and 02:54 UT on December 
13, 2006). The obtained radiative energy during the flare time is 9.04xl0 30 erg. 

It has already been shown that the ejected CME on Dec 13, 2006 carried a me- 
chanical energy of (1.4-4.5) x 10 32 erg. The radiative energy loss in the flare plasma is 
9.04xl0 30 erg. Adding the two shows that (1.49-4.59) x 10 30 erg of energy is utilized. In 
this energy estimation we do not consider the energy used in the production of SEPs, 
magnetic field reconfiguration or unerupted plasma transportation. 

4. Discussion 

The magnetic energy can build up in the process of shearing motion of the magnetic 
footpoints or during the flux emergence processes. The emerging flux region is the location 
for many observed solar flares and CMEs (Martin et al. 1982). The emergence of a 
magnetic field can occur on a small scale (Berlicki et al. 2004) or large scale (Li et al. 
2000). In the case of NOAA AR 10930, a large scale magnetic field emerged near the pre 
existing large south (negative) polarity region. From the magnetic map and transverse 
vectors it appears that the region between the two polarities is highly sheared (Schrijver 
et al. 2008). The flux tubes emerge with a twist (Longcope and Wclsch 2000) albeit with 
various quantities. While the twisted flux tube is emerging, it imports magnetic energy 
and helicity into the corona (Magara and Longcope 2003). In the emerging flux region 
the energy contribution comes from both the horizontal motion of the footpoint as well 
as from the vertical motion of emergence. In our study, we do not know how much is the 
contribution from each of these components. This will be the subject matter of future 
studies. For the present case, however, our results establish that in an emerging active 
region such as NOAA 10930 discussed here, the magnetic free energy was increasing until 
the onset of X3.4 class flare and it accumulated 5.35 xlO 32 erg of magnetic free energy 
over a period of 4 days. 

There have been several studies on this active region to estimate the energy released 
during the X3.4 class flare on December 13, 2006 based on NLFFF model. Schrijver et al. 
(2008) have found that there is a drop in energy of 3xl0 32 ergs during the X3.4 class 
flare. Guo et al. (2008) have obtained a 2.4xl0 31 ergs of drop in energy after X3.4 class 
flare. On the other hand Jing et al. (2010) did not find any drop in energy, instead they 
observed an increase in free energy after the flare. In the present study, we have used 
a magnetic virial theorem to estimate the magnetic free energy available in the active 
region. Even though the cadence of the vector field measurement is poor, we were able 
to see a decrease in the magnetic free energy that persisted even 10 hrs after the flare. 

We have identified a reduction of 3.11 x 10 32 erg of magnetic free energy in AR 10930. 
The CME associated with the energetic event that took place on December 13, 2006 
carried about (1.4-4.5) x 10 32 erg (using projected and deprojected method) of energy. 
The estimated radiative energy loss during the X3.4 class flare was about lxlO 31 ergs. 
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The range of energies gives rise to two possibilities: (1) the entire energy for the CME 
came from the same active region; or (2) only a portion of the CME energy has been 
supplied by this active region. The latter is based on the assumption that the active 
region is associated with a single footpoint of the CME (Simnett and Harrison 1984, 
1985, Harrison and Simnett 1984, Harrison 1986). If the calculated energy is closer to the 
projected calculations (the lower limit) then the former is likely true, while if the energy 
is closer to the de-projected values then the latter applies. 

The latter possibility is more likely for this event given that it is a halo CME and 
therefore likely to be highly projected in the LASCO images. Let us consider the upper 
extreme of the energy. Half of the CME energy is around 2.25 xlO 32 erg, which leaves 
8.6 xlO 31 erg in the lost magnetic free energy available for other purposes. Emslie et al. 
(2004) determined the energy budget on two flare/CME events on 21 April and 23 July 
2002. Both were associated with an X-class flare (XI. 5 for the former and X4.8 for the 
latter). They found that the flares contributed between 1-4 xlO 31 erg, which included the 
energy for particle acceleration. Given that flares are broadband phenomena the emission 
energy contribution extends across the electromagnetic spectrum, and so a larger energy 
budget for emission is expected. We must also consider the energy in moving around solar 
plasma and un-erupted magnetic fields (e.g. Webb et al. 1980). It is not unreasonable to 
conclude that this energy budget combined could accommodate the remaining magnetic 
energy (<~ 8.6xl0 31 erg) that did not contribute to the CME. 

At this stage, however, we do not have sufficient evidence to determine which pos- 
sibility is correct, but it is noteworthy that the calculated energy range (for the CME) 
includes the value for the available free energy, which lies almost exactly halfway (by 
ratio) between the two extremes. 

The actual magnetic free energy release could also be larger than that computed. 
This is because of the following reasons: 

1. The magnetic flux was still emerging before the flare (Park et al. 2010) and our 
preflare magnetic field measurement was about 5 hrs before the flare. So, the 
available free energy could be larger than estimated here. 

2. The post-flare energy calculation was performed at a time one hour after the post- 
flare time. Within this 1 hr gap the free energy would have increased as the flux 
was still increasing (Park et al. 2010). Hence the measured dip in free energy could 
be larger. 

3. The photosphere is not force-free. To make it close to the force-free we have used the 
pre-processing technique that may reduce the actual value of the available energy. 
Metcalf et al. (2008) have applied a realistic model to test the Lorentz forces in the 
photosphere. In their model they compared the net force at various heights. They 
have found that the net force falls close to zero at the chromospheric height but 
not in the photosphere. Due to this net force at the photosphcric level there is an 
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underestimation of the virial energy by about 10%. Hence, the free energy will also 
be underestimated by the same amount. However, this underestimation will not 
affect significantly the change in free magnetic energy measurement between before 
and after the flare. 

4. It is not well known how much the 180° ambiguity will affect the virial energy 
estimation. 

We have identified the source of uncertainty in the energy estimation using the virial 
theorem. The estimated uncertainty is much smaller than the free energy. Metcalf, 
Leka and Mickey (2005) have estimated the uncertainty in magnetic energy estimation 
by pseudo-Monte Carlo method by displacing the origin to different locations on the 
vector fields. This method is valid only if the measured field is forced. However, the pre- 
processing technique makes these fields close to force-free and hence measuring the mean 
value of the energy by shifting the origin to different locations on the image is not required. 
Wheatland and Metcalf (2006) have introduced a new method to find the magnetic energy 
using the virial theorem. However, that method depends on the choice of the model to 
compute the vertical field gradients. In the absence of chromospheric magnetic field 
measurements they have used the linear force-free parameter (a) to compute the vertical 
field gradients. The distribution of magnetic field deviates from the constant a in most 
of the locations in the active region. Hence, a detailed study is required to find further 
details on how this method is superior to the other methods. Until then the method of 
using the pre-processed data for the virial energy estimation is a better choice than those 
methods. 

The largest source of error in the CME calculations is the 3-D information provided 
for the de- projection. We used the location of the flare to provide values for 0, <fi (distance) 
and x (mass) but this assumption is limited. Firstly, it is not appropriate to assume that 
the CME originated from the flare or even the active region, and secondly the CME 
cannot be regarded as a point source (as assumed for the distance measurements). The 
CME is a large structure and for the larger ones the associated flare and active region are 
typically associated with a single footpoint only. In this case, 3-D assumptions based on 
the flare and AR would provide information on just that footpoint. This is also the reason 
why we assume for the maximum case that only half of the CME's energy was provided 
by the AR. To our knowledge no systematic empirical study of the source of energy for 
the CME has yet been conducted so it is unknown what proportion of its energy budget 
arises from the AR. Given the values we arrived at from the present study, however 
limited, it seems that a 50% contribution may be a reasonable first-order assumption. 
Such knowledge may lead to a greater understanding of the onset mechanism for CMEs. 

Jing ct al. (2009) have observed for many active regions that there is a decrease in 
magnetic free energy 15 min before the peak time of the associated non-thermal flare 
emission. Such studies cannot be taken up with the present data as the cadence of vector 
magnetic field measurement is poor. The Solar Dynamic Observatory (SDO) (launched 
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February 2010) provides vector field measurements at high cadence. Hence, the SDO 
data set may provide a complete picture on how the free energy changes close to the 
flare time and during the flare. However, for a better estimation of the available free 
energy we need to have chromospheric vector magnetic field data (Metcalf, et al. 2005). 
Currently, the chromospheric vector field data are not available on regular basis. Until 
they are available, photospheric vector magnetogram with the pre-processing techniques 
will provide the good boundary data for estimating the magnetic free energy. 

5. Conclusions 

Solar flares and CMEs are responsible for the release of large amounts of energy from 
the Sun. During solar flares energy is redistributed partly in the form of broadband 
emission from X-rays to the radio band. The energy is also carried by the thermal 
plasma, accelerated electrons, etc. The vast majority of the energy budget, however, is 
allocated to mechanical energy in the form of an expanding CME. 

Using the high resolution time sequence of vector magnetograms from SOT/SP we 
studied the temporal evolution of magnetic free energy in an active region NO A A 10930. 
We then compared the available magnetic free energy with the energy carried away by 
the halo CME. From the analysis we arrive at the following conclusions: 

1. The free magnetic energy increases with time and about 3.7xl0 32 erg of magnetic 
free energy is accumulated over a period of about 20 hrs before the X3.4 class flare 
and CME launch occurred. 

2. There is a decrease in magnetic free energy of about 3.11 xlO 32 erg during the 
flarc/CME. 

3. The estimated energy carried away by the CME is (1.4-4.5) x 10 32 ergs (using pro- 
jected and de-projected values). This is 0.5-1.5 times larger than the estimated 
magnetic free energy 

4. The estimated radiative energy loss during the X3.4 class flare is 9.04xl0 30 erg. 

We believe that more energy was carried by the CME than was available in the 
active region, as some of the free energy must be allocated to radiative energy loss, 
particle acceleration and plasma and magnetic field reorientation. Given the range 
of energies calculated for CME and the possibility that the free magnetic energy 
may be larger, one could easily conclude that the entire energy was provided by the 
active region. 
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